Jet Shapes in Opaque Media 
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We present general arguments, based on medium-induced radiative energy loss, which reproduce 
the non-gaussian shapes of away-side di-jet azimuthal correlations found in nucleus-nucleus collisions 
at RHIC. A rather simple generalization of the Sudakov form factors to opaque media allowing an 
effective description of the experimental data is proposed. 
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Jet quenching, the strong modification of the spec- 
trum of particles produced at high transverse momenta in 
nucleus-nucleus collisions, is a solid experimental result 
of the RHIC program |l|. At present its main observ- 
able is the strong suppression of the inclusive particle 
yield in the whole available range of p± < 20 GeV/c 
00. This suppression can be understood, in general 
terms, as due to the energy loss of highly energetic par- 
tons traversing the medium created in the collision. To 
further unravel the dynamical mechanism underlying this 
effect the questions of where and how the energy "is lost" 
are of special relevance. The most promising experimen- 
tal probes to investigate these issues are clearly related 
to the modification of the jet structures In the most 
successful approach to explain jet quenching, the degra- 
dation of the energy of the parent parton produced in a 
elementary hard process, is due to the medium-induced 
radiation of soft gluons and the energy trans- 

fer to the medium (recoil) is neglected. As usual, the 
medium-modified parton-shower will eventually convert 
into a hadron jet. In the opposite limit, one could assume 
that a large fraction of the original parton energy is trans- 
ferred to the medium, with a fast local equilibration, and 
diffused through sound and/or dispersive modes. The 
latter possibility has been advocated [ljj as the origin of 
the striking non-gaussian shape of the azimuthal distribu- 
tions in the opposite direction to the trigger particle . 
In this paper we show that the same perturbative mech- 
anism able to describe the inclusive suppression data, 
namely, the so-called radiative energy loss, can account 
for the experimentally observed two-peak shapes in the 
azimuthal correlations. The central point of our argu- 
ment is the need of a more exclusive treatment of the dis- 
tributions to describe experimental data with restrictive 
kinematical constrains. To this end, we supplement the 
medium-modified jet formation formalism with the Su- 
dakov suppression form factor following the well known 
in-vacuum approach. Furthermore, we assume that the 
final hadronic distributions follow the parton level ones 
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(parton- hadron duality). 

The medium-induced gluon radiation. Let us sup- 
pose that a very energetic parton propagating through 
a medium of length L emits a gluon with energy u> and 
transverse momentum k± with respect to the fast par- 
ton. The formation time of the gluon is t[ olm ~ dr- Its 
typical transverse momentum is 
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q is the transport coefficient [4j which characterizes the 
medium-induced transverse momentum squared (<T^) m ed 
transferred to the projectile per unit path length A. In 
this picture the gluon acquires transverse momentum be- 
cause of the Brownian-motion in the transverse plane 
due to multiple soft scatterings during if or m- The typ- 
ical emission angle 



sin 9 



(2) 



defines a minimum emission energy uj ~ {2qy/ 3 below 
which the radiation is suppressed by formation time ef- 
fects |9(. The latter is a crucial observation for the dis- 
cussion to follow. Notice that for energies smaller than u> 
the angular distribution of the medium-induced emitted 
gluons peaks at large values. In 0] a fit to the inclusive 
high-p^ particle suppression gives the values q ~ 5. ..15 
GeV 2 /fm for the most central AuAu collisions at RHIC, 
i.e., u ~ 3 GeV [we have taken cD ~ 2(2q) 1 ^ 3 as indicated 
by numerical results]. This situation can be computed 
from the spectrum of the medium-induced gluon radia- 
tion in the limit k\ < \/2uJq, ui -C u> c = \qL 2 . This gives 
in the multiple soft scattering approximation 
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Compared to the corresponding spectrum in the vacuum, 
the one in © is softer and the typical emission angles 
larger, producing a broadening of the jet signals. 

The parton shower evolution. Equation gives the 
inclusive spectrum of gluons emitted by a high-energy 
parton traversing a medium. For practical applications, 
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however, more exclusive distributions, giving the proba- 
bility of one, two... emissions are needed. How to con- 
struct such probabilities, using Sudakov form factors, is 
a well known procedure in the vacuum. In the medium, a 
first attempt to deal with the strong trigger bias effects in 
inclusive particle production, was that of using an inde- 
pendent gluon emission approximation with correspond- 
ing Poissonian probabilities 13]. Here we propose to im- 
prove this assumption by including the virtuality evolu- 
tion through medium-modified Sudakov form factors. It 
is encouraging to notice that the original formulation is 
recovered once the virtuality evolution is neglected. 

Since we are interested in angular distributions the par- 
ton shower description proposed in 14] is particularly 
convenient. Let us introduce the evolution variable 
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to define the branching of a particle (gluon in general) 
with virtuality q and energy u> in two particles of virtu- 
ality q%, q2 < q and energies u>x = zu>, u>2 = (1 — z)u>. If 
u>X, ll>2 S> mi, TO2, as we will assume, £ ~ 1 — cos #12, 
with #12 = 9\ + 62- 0i, 02 are the angles formed by the 
daughter partons with the parent parton. Consequently 
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i.e., the virtuality evolution can be converted into an evo- 
lution in the variable £. The corresponding probability 
distribution for one branching is 

dP(& Z) = ^ dz ^-P(z) A(£ max ,£)0(£ max -£)0(e-£min) 

(6) 

where P(z) is the splitting function. The Sudakov form 
factor controlling the evolution is 0] : 

A(£, W ) = exp j- ^jy d z^- P(z) I (7) 

and e = Qo/^V? with Qq a cut-off. A(£, oS) can be in- 
terpreted as the probability of no branching between the 
scales £ and £ ma x. 

We propose to generalize the above formulation to the 
in- medium case by noticing that in our notation 6] 9] 
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We define the corresponding equations ©, for the 
medium by changing dl vac / 1 dzdk\ to dP ncd / 'dzdk\ given 
by Eq. ©. We have clear motivations to make this 
Ansatz: (i) it provides a clear probabilistic interpreta- 
tion with the right limit to the most used "quenching 
weights" when the virtuality is ignored; (ii) the 

evolution equations in the case of nuclear fragmentation 



functions in DIS can be written as usual DGLAP equa- 
tions with medium-modified splitting functions P(z) — > 
P(z) + AP(z) 8]. In our case Eq. would correspond 
to AP(z) with the appropriate factors. 

Results. In what follows, we present some simple an- 
alytical estimates following the approach just sketched. 
The general case of arbitrary oj\ and u>2 is difficult to 
study analytically. Let us study, instead, two extreme 
cases: (1) one of the particles takes most of the incom- 
ing energy u>i ^> L02 and, hence, 0i ~ - we name this 
as "J-configuration" ; (2) the two particles share equally 
the available energy, oji ~ L02 and 0i ~ 2 - named as 
"Y-configuration" . 

Let us first consider the case (1) where £ = 1— cos0i2 = 
1 — cos0 with sin0 = k±/ui. From J3J we get 
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Therefore the corresponding Sudakov form factor reads 
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where the small contribution form the integration in z has 
been neglected [we have checked that the results does not 
vary as long asw > Qq] ; a s = 1/3 will be assumed in the 
following. Taking £ = 1 — cos0, £ max = 1 and inserting 
into the probability of splitting, Eq. ©, we get: 

dV{9,z) a s C R , \ a s C R 2 
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Eq. (fTT|l is the probability distribution for a parent par- 
ton to split just once, emitting a gluon at angle with 
fraction z of the incoming momentum. 

The variables of (|llf> are polar coordinates with re- 
spect to the jet axis (the direction of the parent parton) . 
Assuming that this parent parton was produced at 90° 
in the center of mass frame of the collision, we can trans- 
form the coordinates of the emitted gluon to the labo- 
ratory polar coordinates <E> and 0i a b with respect to the 
beam direction z. Normally, one uses the pseudorapidity 
r] = —log tan (0i a b/2). The jabobian is 
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where (3 was integrated in the previous expressions to 
give 27T in the spectra. Taking, for simplicity -q — 0, in 
the most favorable detection region, the answer is simply 
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FIG. 1: The probability of just one splitting as a function 
of the laboratory azimuthal angle A$ for a gluon jet of Ej e t = 
7 GeV. Different medium lengths are plotted for the J- and 
Y-configurations (solid and dotted lines respectively). 
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FIG. 2: Position of the peaks of the A$-distribution and 
comparison with PHENIX data from Ref. 18] for the J- and 
Y-configurations. We also show the analytical estimate (I14H 
without any smearing. 



giving the probability of one splitting as a function of $. 
Thus we reach our objective: the possibility of describ- 
ing non-trivial angular dependences, as shown in Eas. (|ll|l 
and Ijl 3(1 with a perturbative mechanisms. The distribu- 
tion found has two maxima whose positions are deter- 
mined by: 
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The angular shape found in ((13(1 is very similar to the one 
found experimentally. We do not intend to perform a de- 
tailed calculation of the experimental situation. We just 
improve our calculation by introducing a simple model 
to take into account the additional smearing of the jet 
shape introduced by the triggering conditions. Setting 
rj = for the trigger particle, we take into account (i) 
the uncertainty due to the boosted center of mass frame 
of the partonic collision b y in tegrating the probability 
((13(1 in 2A?y with Ar] = 1 [22j; (ii) an additional uncer- 
tainty in azimuthal angle A<I> given by a gaussian with 
a= 0.4 471] ■ Specifically, we take 
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N = 2Ar\\j2-Ku 2 being a normalization factor. The re- 
sults are plotted in Fig. |2] for three different medium 
lengths and £j c t = 7 GeV [the quoted value q = 10 
GeV 2 /fm assures that the results hold for gluon energies 
lj < oj ~ 3 GeV]. In order to compare with the central- 
ity dependence of the position of the maxima measured 
experimentally we simply take L — N^t- Although 
this geometric procedure is too simplistic, it serves us 
to have a feeling of the description of the data plotted 
in Fig. |21 We must also point out that the centrality 
dependence of the transport coefficient is not taken into 



account in these figures. The reduction of q with cen- 
trality (g ~ dN/dy ~ N plir t) makes the radiation more 
collinear when g 1 / 3 < u>, i.e., the radiation will be more 
collinear than plotted for decreasing centrality at fixed 
lo — p ^ 

Let us now consider the case of the Y-configuration, 
where w\ ~ 0J2 and the two angles are similar B\ ~ #2 = 
6. Now the variable £ = 1 — cos #12 = 1 — cos 28 and 
k±/uj 2 = sin 2 9 = ^/2. Changing variables we get 
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Repeating the same procedure followed before 
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In this case, the maxima are outside the borders of the 
physical phase space but still a minimum at $ = oc- 
cours. Anyway, the smeared distribution does present 
maxima. The corresponding curves and positions of the 
maxima are plotted in Figs. and [21 

The picture which emerges from our analysis is the 
following: the splitting probability of a highly energetic 
parton produced inside the medium by a hard process 
presents well defined maxima in the laboratory azimuthal 
angle when the emitted gluons have energies u) < lj, with 
Cj ~ 3 GeV for the most central collisions at RHIC. This 
is a reflection of the jet broadening predicted long ago 
as a consequence of the medium-induced gluon radiation 
0-0. In the case that the experimental triggering con- 
ditions are such that only a small number of splittings 
are possible (by restrictive kinematical constrains, e.g. 
passoc ^ p^ lgg ) these structures should be observed. Our 
estimate, taking into account the parton-hadron duality, 
is then p^ soc < Cj ~ 3 GeV. When, on the contrary, 
the allowed number of splittings is larger, we expect this 
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angular structure to disappear, filling the dip located in 
correspondence of the jet axis. The reason is that the 
inclusive spectrum (which does not present a two-peak 
structure) should be recovered in these conditions. This 
could explain the absence of a two-peak structure when 
the kinematical cut on p^ soc is relaxed (or/and when 
ptrigg j g i ar g er f or fixed p^ 00 ). In this situation the 
broadening is, however, still present with a flatter distri- 
bution which extends in a large range of A<i>. A realistic 
comparison with experimental data would need a much 
more sophisticated analysis, including not only the prob- 
ability of multiple splittings but also the hadronization 
here ignored. Our results are very encouraging though as 
the distributions we obtain for the single splitting resem- 
ble very much the experimental findings, both in shape 
and in centrality dependence. 

Let us also comment on the relation with other ap- 
proaches. In 0| the large angle medium-induced gluon 
radiation has been considered and modifications of the 
angular distribution computed. The obtained jet struc- 
tures are, however, basically gaussian in azimuthal angle, 
the effect of the medium being a broadening on the typ- 
ical jet width. Two main differences with our approach 
are the source of this discrepancy: First, and most im- 
portantly, the angular distributions in [19| are given by 
the inclusive spectrum, which does not present the two- 
peak structure in A<!> - the main point of the present 
paper is to provide a proposal on how to improve this 
description by the inclusion of Sudakov form factors; Sec- 
ond, in the single hard approximation used in 0] for the 
medium-induced gluon radiation the typical emission an- 
gle decreases with increasing in-medium path-length |20| 
as sin 6 ~ 1/ \f~L, making the radiation more and more 
collinear with increasing centrality (in contrast with (J2J 
which is independent on the centrality). 

On the other hand, one of the most popular expla- 
nations for the non-gaussian shape of the away-side jet- 
signals found at RHIC is in terms of shock waves pro- 
duced by the highly energetic particle into the medium. 
In this picture, a large amount of the energy lost must 
be transferred to the medium which thermalizes almost 
instantaneously. In general the released energy excites 
both sound and dispersive modes, and only the first ones 
produce the desired cone-signal (in the dispersive mode, 
the energy travels basically collinear with the jet). The 
energy deposition needed for the sound modes to become 
visible in the spectrum has been found to be quite large 
[2l|. To our knowledge, no attempt has been made so 
far to describe the centrality dependence of the shape 
of the azimuthal correlations in this approach. Given 
the fact that the two formalisms described above rely 
on completely different hypotheses, finding experimental 
observables which could distinguish between them is cer- 
tainly an issue which deserves further investigation. New 
data on three particle correlations are expected to shed 



some light on the problem. Here, we just notice that our 
Y-configurations would lead to similar signatures as the 
ones from the shock wave model. In the most general 
case, however, different configurations, not considered in 
our simple analysis, and characterized by different radi- 
ation angles for both gluons, would produce a smeared 
signal. 
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